MP0065 



LDPC DECODER AND METHOD THEREOF 

Inventors: Zining Wu 

Gregory Burd 

5 



CROSS REFERENCE TO RELATED APPLICATIONS 

The present invention claims priority under 35 U.S.C. §119 (e) from U.S. 
provisional application serial no. 60/214781, entitled "Address Generator for LDPC 
Encoder and Decoder and Method Thereof," filed June 28, 2000, the contents of which 
are incorporated herein by reference. 

The present invention is a continuation-in-part application of application Serial 
No. 09/559186, entitled "Multi-Mode Iterative Detector" and filed on April 27, 2000, the 
contents of which are incorporated herein by reference, 

The present invention is related to the following commonly-assigned, copending 
applications: 

"LDPC Encoder and Method Thereof, filed on even date and assigned 

application Serial No. (Attorney Docket No. MP0064), the contents of 

which are incorporated herein by reference, 

"Address Generator for LDPC Encoder and Decoder and Method Thereof 

filed on even date and assigned application Serial No. (Attorney 

Docket No. MP0063), the contents of which are incorporated herein by reference, 
and 

"Parity Check Matrix and Method of Forming Thereof, filed on even date 

and assigned application Serial No. (Attorney Docket No. MP0069), 

the contents of which are incorporated herein by reference. 
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BACKGROUND OF THE INVENTION 
Field of the Invention 

The present invention relates generally to a linear block decoder in a data 
transmission system. More particularly, the present invention relates to a low density 
5 parity-check code (LDPC) decoder a read channel in a disk drive system. 

Description of the Related Art 

Fig. 1 illustrates a conventional digital data transmission system. As shown 
^ therein, a digital data transmission system comprises a transmitting section 300 for 
,q transmitting user data to receiver 500 via communication channel 401. 

;jjl0 The operation of transmission section 300 will now be explained. Prior to processing 
i by transmitting section 300, input or user data maybe encoded with an error correcting 
Q code, such as the Reed/Solomon code, or run length limited encoded (RLL) or a 
;* combination thereof by encoder 302. The encoded output encoder 302 is then 
w interleaved by deinterleaver 308 for input to linear block code encoder 304 which 
JL5 generates parity data in a known manner utilizing linear block codes. One example of 
4 a linear block code -is a low-density parity-check code (LDPC) which is discussed by 
Robert G. Gallager in Low-Density Parity-Check Codes, 1963, M.I.T. Press and by 
Zining Wu in Coding and Iterative Detection For Magnetic Recording Channels, 2000, 
Kluwer Academic Publishers, the contents of each of which are incorporated in their 
20 entirety by reference. Deinterleaver 308 permutes the data so that the same data is 
reordered before encoding by linear block code encoder 304. By permuting or 
redistributing the data, interleaver 306 attempts to reduce the number of nearest 
neighbors of small distance error events. User data at the output of encoder 302 is 
referred to as being in the channel domain; that is the order in which data is 
25 transmitted through the channel. The order of data processed by deinterleaver 308 is 
referred to as being in the linear block code domain. The parity data from linear block 
code encoder 304 is combined with the data encoded by encoder 302 by multiplexer 306 
for input to channel transmitter 310. 

Transmitter 310 transmits the combined user and parity data from multiplexer 306 
30 typically as an analog signal over communication channel 401 in the channel domain. 
Communication channel 401 may include any wireless, wire, optical and the like 
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communication medium. Receiver 500 comprises a front-end circuit 502 comprising 
analog to digital and equalization circuits. The digital signal front-end circuit 502 is 
input to soft channel decoder 504, which provides probability information of the 
detected data. Soft channel decoder may be implemented by a Soft Viterbi Detector or 
the like. The output of the soft channel decoder 504, which is in the channel domain, is 
converted into the linear block code domain by deinterleaver 510. Deinterleaver 510 is 
constructed similarly to deinterleaver 308. Soft linear block code decoder 506 utilizes 
this information and the parity bits to decode the received data. One output of soft 
linear block code decoder 506 is fed back to soft channel decoder 504 via interleaver 
512, which converts data in the linear block code domain to the channel domain. 
Deinterleaver 512 is constructed to perform the reverse operations of deinterleaver 
510. Soft channel decoder 504 and soft linear block code decoder 506 operate in an 
iterative manner to decode the detected data. 

The other output of soft linear block code decoder 506' is converted from the 
linear block domain to the channel domain by interleaver 514. Deinterleaver 514 is 
constructed similarly to interleaver 512. The output of interleaver 514 is passed on for 
further processing to decoder 508. Decoder 508 is implemented to perform the reverse 
operations of encoder 302. 

Summary of the Invention 

According to a first aspect of the present invention, a method of decoding low- 
density parity-check codes is provided comprising the steps of (a) calculating LLrR m |, for 
each parity check equation, at iteration M, in response to step c, (b) decision aided 
equalizing, at iteration i, in response to step a; and (c) calculating LLrQtm, for each 
parity check equation, at iteration i in response to step b. LLrQtm represents 
information from bit node 1 to equation node m, one for each connection, and LLrR m f 
represents information from equation node m to bit node 1, one for each connection. 

According to a second aspect of the present invention, a method of decoding low- 
density parity-check codes is provided further comprising the steps of (d) updating, for 
each parity check equation, the smallest LLrQtm calculated in step c, at iteration i, (e) 
updating, for each parity check equation, the second smallest LLrQtm calculated in step 
c, at iteration i; and (f) updating, for each parity check equation, an overall sign of 
LLrQtm, calculated in step c, at iteration I Step a is responsive to step d, step e and 
step f, and step c is responsive to step a. 
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According to a third aspect of the present invention, step a is calculated as 



follows: 



llrR 



i-l 
ml 
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;J Step c is calculated as follows: 

U j m* * m 

According to a fourth aspect of the present invention, step c is initialized with 
ry soft channel information. 

:; = : According to a fifth aspect of the present invention, the method further 

IK comprises the steps of (g) converting j LlrAPPk | information into hard information b c ,k, 
(h) calculating an equation syndrome for each parity check equation, Si. (i) summing 
each of equation syndrome in step h, (j) outputting the hard information if the sum in 
step i is equal to zero; (k) repeating steps a-c if the sum in step i is not equal to zero, (1) 
determining if i is less than a first predetermined value; (m) repeating steps a-c, if in 
15 step 1, i is less than the first predetermined value; (n) determining if the sum in step i 
is less than a second predetermined value; and (o) outputting the hard information if i 
is at least the first predetermined and the sum in step i is less than the second 
predetermined value. 

According to a sixth aspect of the present invention, the method further 
20 comprises the steps of (p) for each data k, determining if | LlrAPPk | is less than a third 
predetermined value; (q) for each data k, outputting hard information b c ,k if j LlrAPPk | 
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is at least the third predetermined value; and (r) for each data \ outputting soft 
channel information b s ,kif | UrAPPk j is less than the third predetermined value. 

According to a seventh aspect of the present invention, the method further 
comprises the steps of (s) for each data k, determining if a corresponding parity check 
5 equation is violated, (t) for each data k, outputting hard information b c ,k, if the 
corresponding parity check equation is not violated, (u) for each data k, determining if 
| UrAPPk | is less than a third predetermined value, (v) for each data k, outputting hard 
information b Cj k if ) UrAPPk j is at least the third predetermined value, and (w) for each 
data k, outputting soft channel information b s ,k if | UrAPPk j is less than the third 
10 predetermined value. 

*i According to an eighth aspect of the present invention, a computer program for 

y decoding low-density parity-check codes is provided comprising the steps of (a) 

calculating UrRmi, for each parity check equation, at iteration in response to step c, 
fS (b) decision aided equalizing, at iteration i, in response to step a; and (c) calculating 
IS UrQim, for each parity check equation, at iteration i in response to step b. UrQim 

represents information from bit node 1 to equation node m, one for each connection, 
H and UrR m i represents information from equation node m to bit node 1, one for each 
]s connection. 

j According to a ninth aspect of the present invention, a computer program for 

205 decoding low-density parity-check codes is provided further comprising the steps of (d) 
updating, for each parity check equation, the smallest UrQim calculated in step c, at 
iteration i, (e) updating, for each parity check equation, the second smallest UrQim 
calculated in step c, at iteration i; and (f) updating, for each parity check equation, an 
overall sign of UrQim, calculated in step c, at iteration L Step a is responsive to step d, 
25 step e and step f, and step c is responsive to step a. 

According to a tenth aspect of the present invention, step a is calculated as 
follows: 
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Step c is calculated as follows: 



According to an eleventh aspect of the present invention, step c is initialized 
with soft channel information. 

According to a twelfth aspect of the present invention, the computer program 
further comprises the steps of (g) converting |LLrAPP k | information into hard 
information b^, (h) calculating an equation syndrome for each parity check equation, 
si. (i) summing each of equation syndrome in step h, (j) outputting the hard 
information if the sum in step i is equal to zero; (k) repeating steps a-c if the sum in 
step i is not equal to zero, (1) determining if i is less than a first predetermined value; 
(m) repeating steps a-c, if in step 1, i is less than the first predetermined value; (n) 
determining if the sum in step i is less than a second predetermined value; and (o) 
outputting the hard information if i is at least the first predetermined and the sum in 
step i is less than the second predetermined value. 

According to a thirteenth aspect of the present invention, the computer program 
further comprises the steps of (p) for each data k, determining if | llrAPPk | is less than 
a third predetermined value; (q) for each data k, outputting hard information b c ,k if 
| llrAPPk | is at least the third predetermined value; and (r) for each data k, outputting 
soft channel information b s , k if | llrAPPk | is less than the third predetermined value. 

According to a fourteenth aspect of the present invention, the computer program 
further comprises the steps of (s) for each data k, determining if a corresponding parity 
check equation is violated, (t) for each data k, outputting hard information b C( k, if the 
corresponding parity check equation is not violated, (u) for each data k, determining if 
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| LlrAPPk | is less than a third predetermined value, (v) for each data k, outputting hard 
information b c ,k if | LlrAPPk | is at least the third predetermined value, and (w) for each 
data k, outputting soft channel information b S) k if | LlrAPPk | is less than the third 
predetermined value. 

According to a fifteenth aspect of the present invention, a decoder for decoding 
low-density parity-check codes is provided comprising first calculating means for 
calculating LlrRmi, for each parity check equation, at iteration decision aided 
equalizing means for equalizing LlrRmi, at iteration i, in response to the first 
calculating means, and second calculating means for calculating LLrQim, for each parity 
check equation, at iteration i in response to the decision aided equalizing means. 
LLrQim represents information from bit node 1 to equation node m, one for each 
connection, and LlrRmi represents information from equation node m to bit node 1, one 
for each connection, and the first calculating means is responsive to the second 
calculating means. 

According to a sixteenth aspect of the present invention, The decoder, further 
comprises memory means for storing for each parity check equation, the smallest 
LLrQim calculated by the second calculating means, at iteration i, the second smallest 
LLrQim calculated by the second calculating means, at iteration i; and an overall sign of 
LLrQim, calculated by the second calculating means, at iteration i. The first calculating 
means is responsive to the memory means, and the second calculating means is 
responsive to the first calculating means. 

According to a seventeenth aspect of the present invention, the first calculating 
means calculates as follows: 
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The second calculating means calculates as follows: 
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According to an eighteenth aspect of the present invention, the second 
calculating means is initialized with soft channel information. 

According to a nineteenth aspect of the present invention, the decoder, further 
comprises an address generating means for providing indices of each of the parity 
check equations. 

According to a twentieth aspect of the present invention, the second calculating 
means provides an index of the smallest LlrQim. 

According to a twenty-first aspect of the present invention, multiplying means is 
provided for scaling ILrRmi, calculated by the first calculating means. 

According to a twenty-second aspect of the present invention, the decoder 
further comprising slicing means for converting | LLrAPPk | information into hard 
information bc,k ; equation vector means for calculating an equation syndrome for each 
parity check equation Si, summing means for summing each equation syndrome 
calculated by the equation vector means. The hard information is output if the sum by 
the summing means is equal to zero, and the calculations by the first and second 
calculating means are repeated if the sum summed by the summing means is not equal 
to zero. 

According to a twenty-third aspect of the present invention, the decoder further 
comprises first threshold means for determining if i is less than a first predetermined 
value, and second threshold means for determining if the sum by the summing means 
is less than a second predetermined value. The calculations by the first and second 
calculating means are repeated if i is less than the first predetermined value as 
determined by the first threshold means, and the hard information is output if i is at 
least the first predetermined and the sum by the summing means is less than the 
second predetermined value. 

According to a twenty-fourth aspect of the present invention, the decoder further 
comprises third threshold means for determining for each data k if | LLrAPPk | is less 
than a third predetermined value. For each data k, hard information b c ,k is output if 
| LLrAPPk | is at least the third predetermined value as determined by the third 
threshold means; and for each data k, soft channel information b s ,k is output if 
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| LLrAPPk [ is less than the third predetermined value as determined by the third 
threshold means. 

According to a twenty-fifth aspect of the present invention, the decoder further 
comprises judging means for determining for each data k if a corresponding parity 
5 check equation is violated and third threshold means for determining for each data k if 
| LLrAPPk | is less than a third predetermined value. For each data k, hard information 
b c ,k ? is output if the corresponding parity check equation is not violated as determined 
by the judging means, for each data k, hard information b c ,k is output if | LLrAPPk | is at 
least the third predetermined value as determined by the third threshold means; and 
1A wherein for each data k, soft channel information b s ,kis output if | LLrAPPk | is less than 
-Q the third predetermined value as determined by the third threshold means. 

;y According to a twenty-sixth aspect of the present invention, a decoder for 

}^ decoding low-density parity-check codes is provided comprising a first calculator to 
: : j calculate LLrRmi, for each parity check equation, at iteration i-l, a decision aided 
equalizer to equalize LLrRmi, at iteration i, in response to the first calculator; and a 
second calculator to calculate LLrQim ? for each parity check equation, at iteration i in 
*U response to the decision aided equalizer. LLrQim represents information from bit node 1 
to equation node m, one for each connection, and LLrRmi represents information from 
□ equation node m to bit node 1, one for each connection. The first calculator is 
2^ responsive to the second calculator. 

According to a twenty-seventh aspect of the present invention, the decoder 
further comprises a memory to store for each parity check equation, the smallest 
LLrQim calculated by the second calculator, at iteration i; the second smallest LLrQim 
calculated by the second calculator, at iteration i; and an overall sign of LLrQim, 
25 calculated by the second calculator, at iteration i. The first calculator is responsive to 
the memory, and the second calculator is responsive to the first calculator. 

According to a twenty-eighth aspect of the present invention, the first calculator 
calculates as follows: 
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The second calculator calculates as follows: 



f ttt* t m 

m According to a twenty-ninth aspect of the present invention, the second 

51 calculator is initialized with soft channel information. 

% According to a thirtieth aspect of the present invention, the decoder further 

q comprises an address generator for providing indices of each of the parity check 

■ M equations. 

H According to a thirty-first aspect of the present invention, the second calculator 

lft provides an index of the smallest LLrQLm. 

f 2 According to a thirty-second aspect of the present invention, the decoder further 

C3 comprises a multiplier to scale ILrRmi, calculated by the first calculator. 

According to a thirty-third aspect of the present invention, the decoder further 
comprises a sheer to convert | LlrAPPk | information into hard information b c ,k ; , and an 
15 equation vector circuit to calculate an equation syndrome for each parity check 
equation, si. 

According to a thirty-fourth aspect of the present invention, the decoder further 
comprises a summer to sum each equation syndrome calculated by the equation vector 
circuit. The hard information is output if the sum summed by the summer i is equal to 
20 zero, and the calculations by the first and second calculator are repeated if the sum 
summed by the summer is not equal to zero. 

According to a thirty-fifth aspect of the present invention, the decoder further 
comprises a first threshold detector to determine if i is less than a first predetermined 
value; and a second threshold detector to determine if the sum by the summer is less 
25 than a second predetermined value. The calculations by the first and second calculator 
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are repeated if i is less than the first predetermined value as determined by the first 
threshold detector; and the hard information is output if i is at least the first 
predetermined and the sum by the summing means is less than the second 
predetermined value. 

5 According to a thirty-sixth aspect of the present invention, the decoder further 

comprises a third threshold detector to determine for each data k if | llrAPPk | is less 
than a third predetermined value. For each data k, hard information b c ,k is output if 
j llrAPPk | is at least the third predetermined value as determined by the third 
threshold detector; and for each data k, soft channel information b s? k is output if 
10 | llrAPPk | is less than the third predetermined value as determined by the third 
^ threshold detector. 

According to a thirty-seventh aspect of the present invention, the decoder 
H further comprises a judging circuit to determine for each data k if a corresponding 
parity check equation is violated, and a third threshold detector to determine for each 
1$^ data k if | LlrAPPk | is less than a third predetermined value. For each data k, hard 
s information b c ,k is output, if the corresponding parity check equation is not violated as 
^ determined by the judging circuit, and for each data k, outputting hard information b c ,k 
O if j llrAPPk | is at least the third predetermined value as determined by the third 
Jz threshold detector; and for each data k 9 soft channel information b s ,k is output if 
2Qh (llrAPPk) is less than the third predetermined value as determined by the -third 
threshold detector. 

Other objects and attainments together with a fuller understanding of the 
invention will become apparent and appreciated by referring to the following 
description and claims taken in conjunction with the accompanying drawings. 

25 Brief Description of the Drawings 

In the drawings wherein like reference symbols refer to like parts. 

Fig. 1 is a block diagram of a conventional data transmission system; 

Fig. 2 is a block diagram of a data transmission system in accordance with the 
present invention; 

30 Fig. 3 is an example of a parity check matrix in accordance with the present 

invention; 
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Fig. 4 is a block diagram of a memory arrangement of the decoder in accordance 
with the present invention; 

Fig. 5 is a block diagram of the decoder in accordance with the present 
invention; 

Fig. 6 is an example of a factor graph; 

Fig. 7 is an example of a factor graph illustrating the calculation of LlrQim; 

Fig. 8 is an example of a factor graph illustrating the calculation of UrRmi; 

Fig. 9 is a graph of a function for calculating a soft decision; 

Fig. 10 is a block diagram of a decision circuit of the decoder of Fig. 5; 

Fig. 11 is a block diagram of an alternate decision circuit of the decoder of Fig. 5; 

Fig. 12 is a factor graph of the decoding procedure in accordance with the 
present invention; 

Fig. 13 is a flow diagram of the decoding procedure in accordance with the 
present invention; 

Fig. 14 is a flow chart of the decision procedure of the decision circuit of Fig. 10; 

Fig. 15 is a flow chart of the decision procedure of the alternate decision circuit 
of Fig. 11; 

Fig. 16 is a block diagram of a conventional decoding implementation; and 

Fig. 17 is a block diagram of a decoding implementation in accordance with the 
present invention. 

Description of the Preferred Embodiments 

Reference is now made to Fig. 2, which is a block diagram of a data transmission 
system in accordance with the present invention. In general as shown therein, a 
digital data transmission system comprises a transmitting section 300' for 
transmitting user data to receiver 500' via communication channel 401. The inventors 
have observed that a linear block code encoder is not dependent on a position of a bit 
interleaved. Rather the linear block code encoder only requires a list of equations for a 
given bit. In other words, there is no need to process the data in the order defined by 
the interleaves instead data may be processed in the same order as it is written to the 
channel. This can be accomplished by incorporating an address generator to provide 
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an address of the appropriate equation of the linear block code encoder. This principle 
can be similarly applied to the soft linear block decoder. As a result, deinterleaver 308 
of the conventional system is now replaced by address generator 328, and 
deinterleaver 510 is now replaced by address generator 530. Accordingly, there is no 
requirement for the physical interleaving of data in the receiver 500', since the data 
remains in the same order as the order of bits of data in the channel throughout this 
system. The order of bits of data transmitted through the channel is referred to as the 
channel domain. 

The operation of transmission section 300' will now be explained. Prior to 
processing by transmitting section 300', as in the conventional system, input or user 
data maybe encoded with an error correcting code, such as the Reed/Solomon code, or 
run length limited encoded (RLL) or a combination thereof by encoder 302. Addresses 
for the parity equations of linear block code encoder 304 are generated by address 
generator 328 in accordance with an index of the bits of data, the index being 
determined by address generator 328. Address generator 328 is responsive to counter 
730 under the control of controller 740. Controller 740 synchronizes counter 730 to the 
output of encoder 302 so that counter 730 can provide a count of the number of bits in a 
codeword output by encoder 302 and a count of the number of codewords. In the 
preferred embodiment the data block size is 5000 bits. 

Referring again to Fig. 2, linear block code encoder 304 utilizes the user data 
and address from address generator 328 to provide the parity bits to multiplexer 306. 
Address generator 328 is described in more detail in commonly assigned, copending 
patent application entitled "Address Generator for LDPC Encoder and Decoder and 

Method Thereof filed on even date and assigned application Serial No. 

(Attorney Docket No. MP0063), the contents of which are incorporated herein by 
reference. Linear block code encoder 304 is preferably implemented as a low-density 
parity-check code (LDPC) encoder as described in commonly assigned, copending 
patent application entitled "LDPC Encoder and Method Thereof, filed on even date 

and assigned application Serial No. (Attorney Docket No. MP0064), the 

contents of which are incorporated herein by reference. The parity data from linear 
block code encoder 304 is combined with the data encoded by encoder 302 by 
multiplexer 306 for input to channel transmitter 310. In the preferred embodiment, the 
combined data consists of series of a pair parity bits followed by 34 bits of user data. 
This constraint is established by RLL constraint encoder 302. 
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Transmitter 310 transmits the combined user and parity data from multiplexer 
306 typically as an analog signal over communication channel 401 in the channel 
domain. Communication channel 401 may include any wireless, wire, optical, 
magnetic and the like. 

Receiver 500' comprises an analog to front-end circuit 502 comprising analog to 
digital and equalization circuits. The digital signal from front-end circuit 502 is input 
to soft channel decoder 504, which provides soft or probabilistic information of the 
detected data to soft linear block decoder 506. Soft channel decoder may be 
implemented as a Soft Viterbi Detector or the like, and address generator 530 may be 
constructed similarly as address generator 328 in transmission section 300'. In the 
preferred the Soft Viterbi Detector may be implemented with a Soft-Output Viterbi 
Algorithm which is described in J. Hagenauer and P. Hoeher: "A Viterbi algorithm 
with soft-decision outputs and its applications", Proc. IEEE GLOBECOM '90, Dallas, 
Tex., pp. 47.1.1-47.1.7, November 1989, the contents of which are incorporated by 
reference. 

The soft information output by soft channel decoder 504 remains in the channel 
domain and is decoded by soft linear block code decoder 506, in accordance with the 
address of the parity equations generated by address generator 530. Address generator 
530 is responsive to counter 735 under the control of controller 745. Controller 745 
synchronizes counter 735 to the output of soft channel decoder 504 so that counter 735 
can provide a count of the number of bits in a codeword output by soft channel decoder 
504 and a count of the number of codewords. 

Soft linear block code decoder 506 operates in combination with soft channel 
decoder 504 and address generator 530 in an iterative fashion. Soft linear block code 
decoder is preferably implemented as a low-density parity- check code (LDPC) decoder, 
as described in detail hereinbelow. It is noted that since the soft information from soft 
channel decoder 504 to soft linear block code decoder 506 are both in the channel 
domain, thus as noted above, there is no need for any interleavers or deinterleavers in 
receiver 500'. 

After the iterative process has completed, the output of soft linear block code 
decoder 506 is passed on for further processing to decoder 508. Decoder 508 is 
implemented to perform the reverse operations of encoder 302 or correct for any data 
errors. 
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Prior to discussing the construction and operation of the LPDC decoder, 
reference is now made to Fig. 3 for an explanation of the parity check matrix. Fig. 3 
shows only one example of a parity check matrix. Of course, other parity check 
matrices are contemplated. The preferred matrix is 222 rows (or equations) by 5402 
columns, which comprises 220 linearly independent rows (where 5402 =73*74). The 
matrix can be divided into three tiers of equations having 73, 74 and 75 equations, 
respectively. The set of independent rows can be obtained by canceling the last row of 
the second tier and third tier, namely the 147 th row and the 222 nd row. As shown in 
Fig. 3, the following table shows the values of the elements in the matrix: 



Tier 


i** 1 position 


i th position 


1 


lifr = i(mod73) 


0ifr*i(mod73) 


2 


lifr = i(mod74) 


0 if r* i(mod74) 


3 


1 if r = i(mod75) 


0ifr*i(mod75) 



A matrix having 5402 columns can process a maximum LDPC codeword of 5402 
bits. Of course, as will be appreciated by one of ordinary skill in the art, the matrix 
may be truncated to accommodate a smaller block, however the matrix must be at least 
222 x 4366 which is dependent on the constraint of encoder 302. This constraint is for 
example a RLL constraint. The preferred matrix contains no cycles, since a matrix 
having cycles has degraded performance that degrades significantly. With the first 
tier only, the parity check matrix has a D min =2; by adding the second tier, the parity 
check matrix has a Dmin=4; and by adding the third tier, the parity check matrix has a 
Dmin=6. A further description of the parity check matrix is provided in commonly 
assigned, copending application entitled, "Parity Check Matrix and Method of 

Designing Thereof, filed on even date and assigned application Serial No. 

(Attorney Docket No. MP0069), the contents of which are incorporated herein by 
reference. 

The LDPC decoder is preferably implemented using the sum-product algorithm, 
which is described by Zining Wu in Coding and Iterative Detection For Magnetic 
Recording Channels, 2000, Kluwer Academic Publishers, the contents of each of which 
are incorporated in its entirety by reference. 

A linear code is a set of codewords, x, satisfying the matrix equation (1) 
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H«x=0 (1), 
where H is an M x L matrix, and x is a 1 x L vector. 

The parity check matrix for an LDPC encoder/decoder is sparse, that is a small 
portion of the elements being one, all other elements being zero. An example of a 
parity check matrix is shown in equation 2 that corresponds to the factor graph shown 
in Fig. 6. 



H - 



1011 
0111 



(2) 



which defines the following parity-check equations for the codeword x 



xl+x3+x4=0 
x2+x3+x4=0 

Fig. 6 is a factor graph of the parity matrix of equation 2. The factor graph 
contains two types of nodes the bit node {e.g. bl, b2, b3, and b4) and the check nodes 
(e.g. el, e2). Each bit node corresponds to a bit in the codeword, and each check node 
represents a parity-check equation {i.e., a row in the parity check matrix H). Hence, 
the factor graph for an LDPC code with an M x L parity check matrix H contains M 
check nodes and L bit nodes. An edge between a check node and a bit node exists if 
and only if the bit participates in the parity-check equation represented by the check 
node. The factor graph shown in Fig. 3 is "bipartite" in which the nodes can be 
separated to two groups, namely check nodes and bit nodes. Connections are allowed 
only between nodes in different groups. 

A cycle in a factor graph refers to a finite set of connected edges that start and 
end at the same node. The bold lines in Fig. 6 represent a cycle length of four. As can 
be appreciated by one of ordinary skill in the art, four is the shortest cycle length a 
parity check matrix can have. 

In the sum-product decoding algorithm, there are four types of soft information 
circulating in the decoding process: 

1. LLrPt: The soft information from channel APP (a posteriori probability) 
decoder or soft channel decoder 504 for the first iteration or from decision aided 
equalization circuit 856 for subsequent iterations, one for each bit. 
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2. llrQim: The information from bit node 1 to equation node m, one for each 
connection. 

3. llrRmi: The information from equation node m to bit node 1 ? one for each 
connection. 

5 4. LLrAPPt: the overall soft information after each iteration, one for each bit. 

Each iteration performs three steps of operations (superscripts denote iteration 
number): 

' 1. Each bit calculates the information LLrQ that it passes to the 
connecting equations, which is the sum of the extrinsic information LlrR from 
1Q the previous iteration with the channel information LlrP from the current 

J; iteration, excluding the extrinsic information from the same equation, i.e.: 

I llrCW = LlrPiJ + 2 LLrRm'i'" 1 (3) 

f:] m'=£m 

2 Each check equation (each row) calculates the extrinsic information 

E> for each involved bit. The "LLRXOR" (2 ©) denotes the LLR operation discussed 

; r :~ below. 



S llrRmi' = 2 © UrQt'mi (4) 

20 3. Each bit calculates the overall soft information by summing up all 

the UrRmi's and the ILrPi 

LLrAPPJ: = LLrPi' +2 LlrRmi 1 (5) 
m 

The "LLRXOR" operation is defined as: 
25 LLRXOR(y,z) = log(ey+ez)-log(ey +z +e°) (6) 

max(y,z)+log(l+e- 1 y- z I )-max(x+y,0)-log(l+e- 1 y +z I ) 

where log (1+e-H) can be implemented with a look-up table. Such an operation 
requires 5 additions, 2 absolute values and two table look-ups from a look-up table. 
30 This procedure is referred to as the "Sum Product" (SP) algorithm. 
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The "LLRXOR" operation can be further simplified by eliminating the table look- 
ups, i.e. omitting the log (l+e-H) term in equation (6). The equation is simplified as 
follows: 

LLRXOR(y,z)=-sign(y)sign(z)min( |y|,|z|) (7). 
5 As is apparent to one of ordinary skill in the art, equation (7) requires only 

absolute value and comparison operations. This operation is referred to as the "Signed 
Sum Product" (SSP) algorithm. 

In the SSP algorithm, the magnitudes of the extrinsic information conveyed in 
each equation are determined by the two log-likelihood ratios with the smallest 
10 amplitudes. Other LLR's only contribute to the sign of the information, not to the 
iz amplitude. Fig. 7 shows the llrQim's passed from bits to an equation, and Fig. 8 shows 
the LlrRmi's relayed by the equation. 

:S A conventional implementation of SP or SSP requires L units of memory to store 

gi . the LlrPi:'s; L units of memory to store the LlrAPPi's; ML units of memory to store the 
UrQim's; and ML units of memory to store the LlrRmi's. Soft information is updated bit 
» by bit within each iteration. These implementations are referred to as the "bit-based" 
f .4 implementations. As will be appreciated by one of ordinary skill in the art, the bit- 
;S based implementations require a large amount of memory and increased hardware 
complexity. Fig. 12 illustrates an implementation of SSP, which is generally 
2jfl performed in a serial manner. As shown therein, soft channel decoding is performed 
for the first iteration. Thereafter LlrQim and LlrRm! are sequentially calculated. For 
subsequent iterations i, decision aided equalization is then performed on the 
information. After the decision aided equalization iteration LlrQim is calculated in 
accordance with equation (3), and then UrRmi. is calculated in accordance with 
25 equation (4). The conventional process requires, for a 4000 bit codeword, 
approximately 4000 cycles for soft channel decoding, approximately 4000 cycles for 
calculating LlrQim, approximately 4000 cycles for calculating UrRmi, and approximately 
4000 cycles for performing decision aided equalization. The amount of time required 
by conventional methods degrades performance. 

30 In accordance with the present invention, a pipelined architecture is proposed 

for the SSP circuit, as shown in Fig. 17. In accordance with the present invention, soft 
channel decoding is performed to initialize the SSP circuit, and for iteration i, llrCW is 
calculated. At the same time LLrRmf 1 " 1 for iteration i-1 is also calculated. Based on this 
result and LlrAPPi decision aided equalization is then performed. As can be 

35 appreciated by one of ordinary skill in the art, the pipelined architecture allows for 
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reduced pendency, resulting in higher throughput In accordance with the present 
invention, each iteration takes about 4000 cycles for a 4000 bit codeword. According to 
another aspect of the present invention, an "equation-based" arrangement is proposed 
which implements the same algorithm with much less hardware complexity as 
5 described herein below. 

The "equation-based" arrangement takes advantage of both the sparseness of 
the parity-check matrix and the simplicity in the SSP decoding. It is noted that the 
number of equations for each bit involved corresponds to the column weight of the 
parity check matrix, which in the preferred embodiment is three. Of course other 
1Q appropriate column weights may be employed. According to equation (7), each 
„5 equation node only passes two values to the connected bit nodes if signs of the soft 
information are not considered. Therefore, there is only a need to store two LLR 
:=t magnitudes for each equation, and the signs for all the UrQtm's. for each equation, the 
two smallest LLR's (mini and min2), the position of mini (i), and the sign of the 
t$ equation(s); for each bit are stored. Additionally, three signs are stored corresponding 
to the signs of three UrQim's from the previously iteration. UrRim's can be calculated 
from the stored information as follows: 



UrR 1 ' 

ml 



\ltr(£ 1 Imj 1 mini* ~~^ m , otherwise 



s* • sign 



20 where s m i= (-l) I~I -sign(LLrQLm') is the sign of the mth equation 

1 

(8) 

Referring to block diagram in Fig. 5, the corresponding flow chart in Fig. 13, 
and the corresponding factor graph in Fig. 12, the "equation-based" SSP circuit 
25 performs the following process: 

1. The equation-based SSP circuit is initialized with the output of soft 
channel decoder 504 for the calculation of LlrPi, which is selected by multiplexer 854. 
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During subsequent iterations, multiplexer 854 selects the output of decision aided 
equalization circuit 856, as discussed hereinbelow 

2. At each time clock, address generator 804 calculates the indices of each 
equation for the current bit, and position index generator 802 calculates the position 

5 index thereof. Three LlrRmi's are calculated using Equation 8 by a calculating means or 
calculator or equation 1 LLR update circuit 816, equation 2 LLR update circuit 818, 
and equation 3 LLR update circuit 820. 

3. Each bit calculates the information that it passes to the connecting 
equations by calculating means or calculator as follows, 



p The ILrRmi's are multiplied by a scaling factor by multipliers 836, 838 and 840 

respectively to reduce circulation of any errors. The LLrRmi's are also delayed by delay 
} h circuits 832, 834, and 842, respectively, to synchronize with llrPi. 

T\ 4. Each bit stores the three signs of the LlrQ's in memory means or memory 

1£| sections 708, 712 and 726 for equations 1, 2 and 3 respectively. 

5. Each bit updates the equation information (mini, min2, i, and s) for the 
three equations involving the current bit by equation 1 update circuit 808, equation 2 
update circuit 810, and equation 3 update circuit 830, respectively, in accordance with 
the following procedure: 



itrQj=llrPl+ £ IlrR m , 



(9) 



20 



m 




min 1 * 



else if 



25 



endif 
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5. Each bit calculates the overall soft information by summing up all the 
UrRmi's, and llrPi (via delay circuit 848) as an input to decision aided equalizer 856. 



5 Log-likelihood ratios are converted to probability information by converter 850. 

Converter 850 implements the function shown in Fig. 9. In the preferred embodiment, 

the function is implemented as a piecewise linear approximation , shown as the dashed 
3 lines in Fig. 9. Decision aided equalizer 856 processes LlrAPPi and samples stored in 

samples memory 864 from an equalizer (not shown) in front end 502, Decision aided 
1Q; equalization is described in more detail by Zining Wu in Coding and Iterative Detection 

For Magnetic Recording Channels, 2000, Kluwer Academic Publishers, the contents of 
r] which are incorporated in its entirety by reference, 

s Decision circuit 900, which is responsive to the output of summer 852, address 

I* generator 530 and soft channel decoder 504 determines when to terminate this 
Ijj^ iterative process, as is described in detail below. 

.*« Fig. 4 is a block diagram of a memory arrangement of the decoder in accordance 

\ 3 with the present invention. As shown therein, memory or memory means 750 has a 
partition for each equation, and each equation has a partition for storing magnitude 
data, index data and sign data for iteration i and for iteration M. As used herein, the 

20 term partition refers either to a . physical memory partition or a logical memory 
partition. Equation 1 update circuit 808, for each clock cycle, writes during iteration i 
to memory partition 806a or memory partition 806b i 9 and reads data previously 
stored from iteration i-1 from memory partition 806b or memory partition 806a. 
Multiplexer 902 selects the memory partition for write operations and multiplexer 904 

25 selects the memory partition for read operations. Equation 1 LLR update circuit 816 
reads data previously stored from iteration i-1 from memory partition 806b or memory 
partition 806a as selected by multiplexer 914. Equation 2 update circuit 810, equation 
3 update circuit 830, equation 2 LLR update circuit 818, and equation 3 LLR update 
circuit 820 operate similarly. 

30 Fig. 10 is a block diagram of decision circuit 900, and Fig. 14 is a flow chart 

thereof. Decision circuit 900 is responsive to the output from soft channel decoder 914, 
LlrAPPi from summer or summing means 852 and address generator or address 
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generator means 530. Decision circuit 900 comprises memory 914 for storing hard 
decision data b s ,k (for k=l:5402) from soft channel decoder 504, sheer or slicing means 
888 for converting | llrAPPt | information into hard information b c ,k ? memory 904 for 
storing b c ,k, a threshold detector 902 to determine if | llrAPPi | from summer 852 is less 

5 then a threshold value (t2). Decision circuit or decision means 900 further comprises 
counter 920 for counting the number of iterations performed by the decoder and 
threshold detector 910 to determine if the number of iterations exceeds threshold value 
t3, which in the preferred embodiment is four. Equation vector circuit equation vector 
means 908 determines the equation syndromes, one for each equation, S m , m=l:222 in 
10 which S m =0 when equation i is satisfied or S m =l when equation m is violated. 

*j Summer 906 sums each of the S m m=l:222. 

The operation of decision circuit 900 will now be explained in conjunction with 
the flow chart of Fig. 14. After an iteration of iterative decoding, the equation 
■P syndrome is calculated by equal vector circuit 908 and summed by summer or 
15 summing means 906. Threshold circuit 912 determines if the sum of all S m =0, then a 
~ signal is sent to multiplexer 854 to output b c ,k from memory 904. If the sum of all S m 
Ll ^0, then a determination is made by threshold 910 whether the number of iterations is 
□ less than threshold value t3. If the number of iterations is less than threshold value t3, 
3 the next iteration is performed. Otherwise threshold circuit 912 determines if the sum 
of all Sm is less than threshold ti. If so, then a signal is sent to multiplexer 854 to 
output b c ,k from memory 904. Alternatively, the next step in the process is to 
determine if | llrAPPi | is less then threshold value Oa), if so, multiplexer 854 selects 
b s ,k to be output. If | llrAPPi. | is not less then threshold value (t2>, multiplexer 854 
selects b c ,k to be output. 

25 Fig. 11 is a block diagram of decision circuit 900' and Fig. 15 is a flow chart 

thereof. Decision circuit 900 ? is responsive to the output from soft channel decoder 914, 
LlrAPPt from summer 852 and address generator 530. Decision circuit 900' comprises 
memory 914 for storing hard decision data b Sf k (for k=l:5402) from soft channel decoder 
504, sheer 888 for converting | LlrAPPi | information into hard information b Cj k, memory 

30 904 for storing b c> k, a threshold detector 902 to determine if | LlrAPPi | from summer 
852 is less then a threshold value (t2). Decision circuit 900' further comprises counter 
920 for counting the number of iterations performed by the decoder and threshold 
detector 910 to determine if the number of iterations exceeds threshold value t3, which 
in the preferred embodiment is four. Equation vector circuit 908 determines the 

35 equation syndromes, one for each equation, S m? i=l:222 in which S m =0 when equation i 
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is satisfied or S m =1 when equation i is violated. Summer 906 sums each of the S m 
m=l:222. 

The operation of decision circuit 900 ? will now be explained in conjunction with 
the flow chart of Fig. 15. After an iteration of iterative decoding, the equation 
5 syndrome is calculated and summed by equal vector circuit 908 and summer 906. 
Threshold circuit 912 determines if the sum of all S m =0, then a signal is sent to 
multiplexer 854 to output b c ,k from memory 904. If the sum of all S m then a 
determination is made by threshold 910 whether the number of iterations is less than 
threshold value t3. If the number of iterations is less than threshold value t3, the next 
10 iteration is performed. Otherwise threshold circuit 912 determines if the sum of all S m 
!□ is less than threshold ti. If so, then a signal is sent to multiplexer 854 to output b C) k 
^ from memory 904. Alternatively, the next step in the process is to determine if bk is 
y involved in an equation that is violated. If bk is not involved in an equation that is 
^ violated multiplexer 854 b c ,k is provided as the output. Otherwise threshold 902 
f | determines if | LLrAPPk | is less then threshold value (t2), if so multiplexer 854 selects 
^ b Sj k to be output. If | LLrAPPk | is not less then threshold value (t2), multiplexer 854 
selects b c ,k to be output. 

S While the invention has been described in conjunction with several specific 

S| embodiments, it is evident to those skilled in the art that many further alternatives, 

2H; modifications and variations will be apparent in light of the foregoing description. 
More specifically, while the present invention is preferably implemented as an 
integrated circuit, it is contemplated that the present invention may also be 
implemented as discrete components or a general-purpose processor operated in 
accordance with program code instructions or computer program or combination 

25 thereof. These program code instructions can be obtain from a medium, such as 
network, local area network, the Internet, or storage devices. Such storage devices 
include, by way of example, magnetic storage devices, optical storage devices, 
electronic storage devices, magneto-optical device and the like. Thus, the invention 
described herein is intended to embrace all such alternatives, modifications, 

30 applications and variations as may fall within the spirit and scope of the appended 
claims. 
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